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Abstract 

We examine the performance of a quantum phase gate implemented 
with cold neutral atoms in microtraps, when anharmonic traps are em- 
ployed and the effects of finite temperature are also taken into account. 
Both the anharmonicity and the temperature are found to pose limita- 
tions to the performance of the quantum gate. We present a quantitative 
analysis of the problem and show that the phase gate has a high quality 
performance for the experimental values that are presently or in the near 
future achievable in the laboratory. 

PACS number(s): 03.67.Lx, 32.80. Pj, 34.90. +q 

1 Introduction 

The implementation of quantum logic gates (16) is a major goal in the 
current research in quantum information. Several schemes have been pro- 
posed in the latest years, based on different physical systems: trapped 
ions pP| or neutral atoms 0, cavity-QED and photons molecules 
quantum dots and Josephson junctions 5 . The aim is to implement a 
fundamental logic quantum gate that works as a constituent block of a 
quantum computer 0. One of such gates is the phase gate, whose truth 
table is 
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Atoms are very good candidates for implementing quantum gates, be- 
cause of the significant experimental achievements realized in recent years. 
The techniques to cool and trap charged and neutral atoms have lead to 
an unprecedented precision in controlling even single atoms. In particular, 
neutral atoms seem to be the most promising systems for quantum infor- 
mation processing, because the dissipative influence of the environment is 
relatively weaker when compared to other physical systems. 

A proposal for implementing a phase gate with cold neutral atoms 
stored in microtraps has been recently put forward by Calarco et al. 1 . 
Two atomic internal states, denoted as | a) i — >| 0) and | b) i — >| 1), are 
used as logical states, and the operations that realize the truth table 



Eq.Q involve the external degrees of freedom. Each atom is placed in a 
microtrap, which can be state-selectively switched off and substituted by 
a larger harmonic potential that allows collisional interaction between two 
atoms. The interaction provides the phase that appears in the truth table 
Eq.Q. For the sake of simplicity, the traps were assumed to be perfectly 
harmonic. 

In the present article we reexamine this proposal for a phase gate and 
in contrast to the work in Ref.[3], we employ, when it is possible and useful, 
the exact analytic expression of the eigenstates of two harmonic oscillators 
with contact interaction [SJ- Besides, with respect to the harmonic term, 
we consider the successive terms of the Taylor series expansion of the 
potential. In particular we focus on the fourth term, which yields the 
lowest-order correction to the dynamics. The effects of the temperature 
are also examined. In realistic experimental situations these feature will 
unavoidably become relevant. 

The article is organized as follows: In Sec.[3]we describe how the quan- 
tum phase gate analyzed here can be implemented. We stress that two 
conditions, (i) full revival of the vibrational state and (if) the acquisition 
of the correct phase shift, are the essential ingredients for a correct per- 
formance of the gate. Then section |2~T1 describes how the phase gate can 
be implemented by using neutral atoms in microtraps and how the two 
conditions mentioned above can be fulfilled. In Sec.[3]we examine the per- 
formance of the gate when the atoms are at zero temperature and oscillate 
in harmonic traps. This situation was already investigated numerically in 
but here we use the exact eigenstates of the problem. The results 
presented here are in total agreement with those shown in 0. In Sec.|I] 
we examine the performance of the gate when the trap is anharmonic. We 
obtain quantitative estimates for the quality of the gate performance. In 
Sec. Q3 we present a heuristic method that, given a certain anharmonic- 
ity, improves the performance by choosing the trap parameters in such 
a way as to optimize the overlap between the initial and the final state. 
In Sec. ED we consider the case when the atoms are at finite temperature 
and we provide a quantitative measurement of the performance using a 
definition of the fidelity. In Sec. |7| we show some connections with our 
anharmonic model and the physical implementation of the gate on atom 
chips. Section 03 contains our conclusions. 

2 Implementing a phase gate with neu- 
tral atoms 

A phase gate with the truth table Eq.Q can be implemented by employing 
two internal atomic states (hyperfine states) as the logic values a and b 
and by making use of the vibrational degrees of freedom of the atoms to 
manipulate the qubits. In order to keep the exposition simple, we assume 
that the atoms are at zero temperature, described by the state 

|*A B (i = 0)> = I Vab(O))® I x> 

= | t/>ab(0))(c q I a A ) + c b | bA>)(c' a | a B ) + c' b | 6 B )) 

= | ^Ab(O)} (c a C a | Ha) | d B ) + C a C b | Oa) | &b) 

+c b c' a I &a) I a B ) + c b c' b | &a) I 6b)), (2) 

where 



2 



j x) = (c Q | o A ) + c b | 6 A ))(<4 I ob) + c' b | 6 B >) 



(3) 



is the general initial internal state of the two atoms, the complex coef- 
ficients c a , c b , c' al c' b satisfy the normalization conditions c a 2 + | c b \ 2 — 1 
and | c' a | 2 + | c' b \ 2 — 1, and | V'ab(O)) is the vibrational state at t = 0. 
The phase gate operation is obtained with a sequence of unitary transfor- 
mations that lead to the final state 



at t = t. A comparison between the two expressions Eqs.Q and J2J 
for the initial and final states shows that two ingredients are essential: 
(i) a sign change must occur only in the last term of Eq.Q and (ii) 
the vibrational state must be disentangled from the internal states at 
the end of the gate operations. These two conditions are fulfilled if the 
motional state, whose wavefunction is ip{xA, XB,t), regains its initial form 
iI>(xa, xb, 0) at some later time r and acquires a phase it only when both 
atoms are in the excited state. 

2.1 Phase gate with two trapped cold atoms 

A natural choice to obtain the recurrence of the initial state ip(xA, xb, 0) 
are atoms oscillating in harmonic traps, where full revivals of wave packets 
are periodically observed for any initial state of a single atom. However, 
the sign change, i.e., the occurrence of a phase tt in the wave packet only if 
both atoms are internally excited, can be achieved only via an interaction 
between the atoms that depends on the internal states. This interaction 
provokes deviations from full revivals of the wave packet and works against 
a correct performance of the phase gate. The aim of our studies is to 
investigate under which conditions the complete wave packet revival is 
at least approximately satisfied, and to evaluate the corresponding gate 
fidelity. 

The system we consider to implement logic quantum gates is an array 
of cold bosonic neutral atoms confined in microtraps. We briefly describe 
how the gate works at temperature T = 0. More details can be found in 



We assume that bosonic rubidium atoms are employed and use typical 
experimental values for the parameters. For the sake of simplicity, we 
focus our attention on only two atoms, under the assumption that we 
can restrict our analysis to a one-dimensional system. For this purpose a 
strong harmonic confining potential, with frequency lu±, in the transverse 
directions y and z can be employed. 

According to Fig. at t < the two atoms are confined in two har- 
monic microtraps of frequency luq, centred at x — — xq and x = xo, re- 
spectively. The distance between the two traps is such that the atoms do 
not interact each other. At time t = 0, the shape of the trapping potential 
changes for the particles in the state |6) into a common harmonic well of 
frequency ui < Wq, centred at x = [dashed line in Fig. Q (b)], whereas 
the potentials for the particles in state \a) remain unchanged [solid line in 
Fig.|I[[b)]. By removing the barrier, particles in state |6) start to oscillate 
and will collide. As a last step, the atoms have to be restored to the initial 
motional state of Fig. 0^a). The whole process of switching potentials is 
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performed through switching the shape of the potential instantaneously 
at times t = and t = r, where r is a multiple of the oscillation period 
in the well of Fig. H^b) (dashed line) . 




b) 




Fi gUT6 li Configuration at times t < and t > r (a), and during the gate operation (b). 
The solid (dashed) curves show the potentials for atoms in internal state \a) (|6)). 



In order to avoid undesired interactions between the two atoms in 
different internal states, the atom in the ground state is shifted in the 
transverse direction. Indeed, this interaction would spoil the performance 
of the quantum gate, as already discussed in [2]. Only when both atoms 
are excited they oscillate in the central trap and interact. 

For convenience of notation, let us define the following Hamiltonians 
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-xb) (9) 

The Hamiltonians Q and @ describe the atoms oscillating in the har- 
monic microtraps with frequency wo, centred in xo and — xo, respectively, 
where the two atoms are placed before the gate operation. The atoms 
remain in these traps when they are in the internal ground state. The 
Hamiltonian (J7J describes one atom freely oscillating in the central trap 
with frequency u). The Hamiltonian © describes an anharmonic central 
trap, derived from the harmonic oscillator Hamiltonian JJj by adding a 
quartic term, where the dimensionless parameter A measures the strength 
of the anharmonicity. The choice of this particular form of the anharmonic 
trap will be justified in Sec.|l] Finally, the Hamiltonian Q describes two 
atoms in the anharmonic central trap interacting via a contact potential 
described by the Dirac delta function S(x\ — xb)- The coupling strength 
depends on the three-dimensional scattering length a s of the two atoms in 
the internal |6) and on the frequency uj±_ \12\ I13| . We recall that the ap- 
proximation 2a s huij_8(xA — xb) is valid only for l± = [h/(Muj±)] 1 / 2 3> a s . 
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The initial state of the two atoms is 

I V'ab(O)} =| ^o A )(c a | a A ) + c b | b A » | <Po B )(c a \ a B ) + c b \ b B )) (10) 

where (pi~ denote the eigenstates of H^ 1 . The gate operation is described 
by a unitary evolution operator U a ,p{t), which depends on the internal 
states a, j3 — a, b of the two atoms and transforms the initial state into 

I Vab(*)) = U a , p (t) | Vab(0)) 

= (e _! " ot | ip~ A )) ( e -^ * | <^+ B }) Ca c' a | o A ) I ob) + 

(e"'" ot | ^oa)) (e^ fl B f | Vob>) Cac'i, | a A ) | 6b) + 

(e-*^ 1 ! ¥>oa)) (e^ ' I ^o + b)) ftci | 6 A ) | a B > + 

(e _i ^ Bt I <^oa) I I 6a) I 6 B ) (11) 

at a later time t. Here the anharmonicity of the central trap has been taken 
into account. The state Ea. l|ll^ is in general no longer a separable state of 
motional and internal degrees of freedom. However, in the present scheme 
the separation between the states of the motional and internal degrees of 
freedom can be realized to a good approximation. Indeed, the state 

| ^ABit)) = (e— ot I ^o A » (e- iW0< | <p+ » (12) 

describes two harmonic oscillators in two well-separated microtraps and 
has therefore full revivals with period T® sc = 2-k/uiq for any initial state. 
The states 

| tf B (i)> = (e-*- 1 I <p^)) (e-^ | p+ )) (13) 

and 

| i>%{t)) = (V^ 4 | ^)) (e- to0< | p+ » (14) 

describe one atom in the microtrap and the other in the central trap. In 
this case, the atom in the microtrap is shifted in the transverse direction 
in order to avoid undesired interaction between the atoms. Therefore, if 
the central trap is harmonic (A = 0), full revivals of the two wave packets 
at different periods occur. 
The state 

| ^ B (t)) ee (e-^AB* | V - A ) I (is) 

is affected by the interaction between the atoms in the wide trap. This 
interaction is necessary in order to yield the sign change for the phase 
gate operation but it also modifies the atomic wave packet. For a good 
performance of the quantum gate the modification must be small. The 
overlap fidelity 

0(tf&,t)= (VabW I^ab(0)) (16) 
between the initial vibrational state | V'ab(O)) and that at a later time 
t > 0, | V'abW) gi ves an estimate of the quality of the gate performance. 
If the revival of the motional state is nearly complete at t — r, it results 
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0« b B ,r)| 2 ~l. (17) 
Under this condition, we can write the motional state at time r as 



I VabOO) - e^ 6 ^ | Vab(O)) (18) 
where <f>bb(T~) denotes the phase of the motional wave function. When 

(f>bb(r) = 7r,37r, 5tt, . . . , (19) 

the phase gate operation is correctly realized. In the next sections we in- 
vestigate when the two conditions of full or nearly full revival (expressed 
by Eq. C2J) of the vibrational state of the excited atoms and the acquis- 
tion of the correct phase (expressed by Eq. dJ) are satisfied. 



3 Ideal case: two cold atoms in harmonic 
traps 

We examine first the performance of the quantum phase gate when all 
traps are harmonic. This ideal situation was already investigated in [2], 
but we examine it again, because here we use the exact solutions of the 
Schrodinger equation for this problem. Indeed, it has been recently found 
that the problem of two interacting atoms in one harmonic trap has an 
exact solution in one, two and three dimensions [Hj- Here we simply 
summarize the results for one dimension. It is useful to define the new 
coordinates X = (ia + xb)/V% and x = (a; a — Xs)/V2. The pseudo- 
particle described by the centre of mass coordinate X is a free harmonic 
oscillator, with eigenstates tp„(X) and energy E n — Tiuj(n + 1/2). The 
pseudo-particle described by the relative coordinate x describes a har- 
monic oscillator that feels the contact potential at the origin x — 0. Its 
eigenstates split into two subsets, depending on the function parity. The 
odd eigenstates are still those of the free harmonic oscillator, since the con- 
tact potential acts only at x = 0, where the odd wave functions vanish. 
The even solutions 
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have energy E v = Tuj(2v + 1/2). Here the real parameters v are solution 
of the transcendental equation 



r-(-iz) 



1 a s lo± 



(21) 



where a x = yJh/(Muj) is the characteristic length associated with the 
central trap and the normalization coefficients 



B v 



r(i/2 - ^)r(-z. 



(22) 



j#(1/2 - v) - 4>(-u)] 

are defined with the help of the logarithmic derivative ip of the gamma 
function. 

With the help of these exact eigenstates we evaluate the fidelity Eq. HI 711 
and the phase shift Eq. (119^ when both atoms are excited. At t < the 
two atoms are in the vibrational ground states of their own microtrap, 
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say, atom A in the left trap, centred at x = — xo, and atom B in the right 
trap, centred at x = xo- At t = the microtraps are switched off, the 
central trap is switched on and both atoms oscillate in the same trap and 
interact. In order to calculate the overlap fidelity Eq. I17L we can pro- 
ceed in two different but equivalent ways. Since the atoms are identical, 
the gate performs correctly even when the two atoms end up in the other 
trap at the end of the operation. We can therefore either symmetrize the 
initial wave function for the two bosonic atoms, or estimate the overlap 
fidelity Eq. Ijl7^ as the sum of the probabilities to find the (distinguish- 
able) atoms in the original traps or with the initial positions interchanged. 
Since in this section we use the wave functions of the center of mass X 
and of the relative coordinate x, the first approach is more convenient. 
The two-atom vibrational state at generic time t > is 

^(X,x,t) = ^c„,„e- ,( " +2 " +1) "VnW^W (23) 

where the odd eigenstates of the relative coordinate x are not included 
since they describe fermions. The coefficients c n ,„ vanish when n is odd, 
otherwise 



y n(uj +uj) ( n /2)! \lu + ujoJ 

(24) 

where 



1 V = / dye-( 1+ ^^U (-is, 1/2; y 2 ) cosh 

J — oo 

In Fig. |21 we show the overlap fidelity, Eq. I|17|l . The revival of the 
vibrational state, occurs with periodicity T osc /2. Indeed, after this period 
the two atoms are in x = ±xo- The revival is nearly complete, as the 
overlap fidelity approaches the value 0.99. In Fig. |2]we also report the 
situation in which atoms in different internal states feel a contact potential 
in order to show what happens. 

In Fig.|5]we show the phase shift (/^(r) due to the interaction. There 
is a fast change of the phase between the two revivals, in correspondence 
to the presence of the two atoms at the bottom of the trap, where the 
interaction occurs. The figure suggests that one can assume that at each 
interaction there is a jump in the phase. Therefore, after a suitable number 
of collisions, the vibrational state acquires the correct phase for the gate 
operation. The exact results shown here confirm the validity and accuracy 
of the results of Ref. 2 obtained numerically. 
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4 Phase gate performance with anhar- 
monic traps 

The operation of the phase gate relies on several simplifying assumptions. 
In the previous section we have assumed that the atoms are at zero tem- 
perature and oscillate in harmonic traps. The experimental conditions 
are necessarily less ideal and the problem of estimating the effect of de- 
viations from the ideal conditions is particularly important. Among the 
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Figure 2: Dynamics during gate operation: projection of the initial state on the state 
evolved without (top) and with interaction (center); projection of the evolved state on the 
corresponding state evolved without interaction (bottom left) and interaction-induced phase 
shift (bottom right). We choose u> = 2tt17. 23kHz and luj_ = 2-7rl50kHz with the initial wells 
having a frequency luq = 2uj and displaced by xq = 410 nm. This reproduces analytically the 
results obtained in f2] with a purely numerical approach. 



causes that can lead to bad performance of the phase gate, we mention 
random noise, caused by fluctuating electromagnetic fields, temperature 
effects and anharmonicity of the trapping potentials. In particular, some 
anharmonicity might more easily appear in the wider central trap and 
seems to be the most important disturbance to be taken into account. 

In this section we study the gate performance when the excited atoms 
oscillate in an anharmonic trap, as described by the Hamiltonians Eqs. (|SJ 
and Indeed, independently of its exact expression, we can expand the 
trapping potential in Taylor series. The first correction to the harmonic 
approximation is a cubic term, that we shall neglect since at first order 
of approximation it does not lead to any correction to the energy and it 
does not affect the atomic motion. 

The next relevant correction to the harmonic trap is a quartic term. It 
is symmetric and it is responsible for an energy shift. We neglect the other 
terms in the Taylor expansion, which yield minor perturbations. The atom 
dynamics in the central trap is then described by the Hamiltonians (|HJ 
and (jnj- 

The atomic wave packet does no longer show full revivals. Moreover, 
the partial revivals of the initial state that still occur are no longer pe- 
riodic, in a strict sense. Therefore, the question arises as to how much 
anharmonicity can be tolerated without destroying the good performance 
of the phase gate. 

4.1 Anharmonic trap: Overlap fidelity for the in- 
ternal states | (Za) | &b) and | bx) | ae) 

In this subsection we examine the performance of the gate when only 
one atom is in the internal excited state. During the gate operation the 
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excited atom (atom A in the left trap, say) oscillates freely in the central 
anharmonic trap, while atom B oscillates as a free harmonic oscillator in 
its microtrap. We need only to focus our attention on the motion of atom 
A, whose vibrational state at time t can be expanded on eigenstates of 
the harmonic oscillators according to 

| Va(*)) = ^(t)e~ i2mi ! <PnA) (t — jT^) (26) 

n 

where tp n denotes the eigenstates of Ho and T oac is the period of oscillation 
of the central trap [dashed line in Fig. EJb)]. The expansion coefficients 
Ck(t) satisfy the differential equations 



Cn(t) =-i|A I v/(n + l)(n + 2)(n + 3)(n + A)c n+i {t) e - iM 

+2(2n + 3) v /(n + l)(r l + 2)c n+2 (t)e- l4 ' rt 

+3 [(n+l) 2 + n 2 ] c n (t) 

+2(2n - l)v , (n-l)nc n _ 2 (t) e ' 4 ' r ' 

+ ^/(n - 3)(n - 2)(n - l)nc„_ 4 We*') , (27) 
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Figure 3: Dynamics during gate operation: projection of the initial state on the state evolved 
with different values of anharmonicity. Trap parameters have the same values as in Fig. |2] 



The overlap fidelity is 
O ba (t) | 2 =| <^ A | e-**** | ^oa) I 2 =l ^ Cfc We-^ fct ^(0) j 2 (28) 

k 

where 

Cfc(0) s {<p kA \ ^- A ) (29) 

The overlap fidelity (1161 has been evaluated numerically and is plotted 
in Fig. |3 on the top for different values of the parameter A. It is evident 
from the figure that the gate tolerates some anharmonicity, with a thresh- 
old value of the order of A ~ 10 -4 . The specular case of atom B excited, 
while atom A is not excited, gives the same results. 
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4.2 Anharmonic trap: Overlap fidelity for the in- 
ternal state | 6a) I &b) 

Here we examine the performance of the gate when both atoms are excited. 
We expand the motional state of the two atoms on the eigenstates of the 
central harmonic trap 

k,l ° SC 

and the expansion coefficients Ck t i(t) statisfy the equations 
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+ (4Z + 6)y/(l + + 2)c k , l+2 (t)e- 4M 

+ y/(l + !)(/ + 2)(Z + 3)(Z + 4)c fc , i+4 (t)e- 8 H (31) 



where 



,4 fetam = -===-=== / d£ e " 2 « Hk(£)Hi(£)H n (£)H m (£) 
V2 k + l + n + m VkU\n\m\ J_ (X1 

(32) 

and denotes the Hermite polynomial of order n. The first term in 
Eq. l-'ilt comes from the contact interaction between the two atoms, 
whereas the other terms are due to anharmonicity. Assuming that the 
atoms are distinguishable, the overlap fidelity is 



O bb (t) = ! (<Poa<Pob I e~i H AB l | cpg A ipg B ) | 2 
+ I {<Poa<Pob I e"* H ^s* | <Po A <Pob) I 

E/,\ -ik2irt -U2irt AB*f n \ 
c k i(t)e e C kl (0) 
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where we have defined 



Cki{t)e e c M (0) | (33) 



c M B *{0) = {(PqaVob \ <PkA<PlB) 

ch A *{0) = (v>qa<Pob I <PkA<piB) (34) 
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We have numerically evaluated the overlap fidelity, which is shown in 
Fig. |21 on the bottom for different values of the parameter A. Also in this 
case we see that an anharmonicity of the order of A ~ 10 -4 or less does 
not prejudicate the performance of the phase gate. 

In conclusion of this section, we note that different choices of the val- 
ues of the parameters lead to very different performance qualities. For 
a fixed value of the anharmonicity parameter A different gate perfor- 
mances are obtained for different values of the other parameters. From 
the first term on the right hand side of Eq. il-'-> I li one sees that the effect 
of the contact interaction on the atom dynamics depends on the value 
of uj±/u}a a (Muj/h) 1 / 2 . If this quantity is larger than w 0.7, it spoils the 
occurrence of full revivals; if it is too small, too many oscillations are 
needed to create the phase <fibb(t) = it. We also note that the frequency 
lu±_ must be large enough to prevent excitations along the transverse di- 
rection (Zj_ 3> a a ). In spite of these limitations, it is possible to find 
reasonable values for these parameters that make a correct performance 
possible, as the data in Fig. [3] (bottom) show. 

5 Optimization of gate performance 

The gate perfomance can be optimized reducing the number of oscillations 
and selecting the trap frequencies w{X) and w_l(A) such that the overlaps 
\0 (tpa/3, t)\ are close to one. In this way we have two effects: better 
performance and faster gate. In Fig. 0]we report the overlaps \0 (tpa/3, t)| 
for two different anharmonic situations. 
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Figure 4: Overlaps |o (ip a i3,T} I in order to optimize the performance of the phase gate. 
In Table we report the values that optimize the performance. 
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Table 1 : Trap frequencies maximizing the fidelity curves in Fig. |1] 

With these values we obtain the results shown in Fig. The crosses 
in the bottom pictures of Fig. are obtained by means of these two 
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assumptions: (i) the particles move against each other, come in contact 
during a certain time interval [ti,tp] and then separate again; and (ii) the 
velocity of each particle and the shape of its wave function do not vary 
during the interaction. It follows that (see Ref.[5j for more details) 



4>bb (Tosc) — 2 



ui±a s 



(35) 



in harmonic oscillator units. Here the velocity v is a positive constant 
value given by 

v = |&<l&±(t)|a#±(t)>|«=ij = 



n + 1 fn . . 

— — c n+ i (tk) - Y g Cn-l (tfe) 



+ 



+ c n (/; fc )c 9 (t fc ) e' 



where 2^ +1 „ and „ are given by 



n + 1 T 4 /n-4 

2 -Ln + l,q "T Y o n— 



[2 n 



die 31 x s H n (x)H q (x), 



(36) 



whereas the coefficients c„(7j) are given by equation l|27|l and tfc = (2fc + 1) T OBC /A 
with an integer. 



o 



o 



fP-5 

O 






0.5 
t/T 

osc 




0.5 
t/T „sc 



Figure 5: Dynamics during gate operation: projection of the initial state on the state evolved 
with interaction (first two rows); projection of the evolved state on the corresponding state 
evolved without interaction (third row) and interaction-induced phase shift (bottom). The 
displacement of the two initial separated wells is the same as in Fig. [5] 



Fig. shows that there is a good agreement between the numerical 
result (solid line) and that given by equation 13 5 ^ . The agreement is not 
perfect, though, simply because the velocity is not constant during the 
interaction. 
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6 Gate performance in anharmonic traps 
at finite temperature 

Now we examine the gate performance when both anharmonicity and tem- 
perature effects are taken into account. Since the atoms are not in the 
ground state of the initial trapping potential and therefore their temper- 
ature T 7^ 0, it follows that, for a finite temperature T, the initial state 
of the two atoms is described by the density matrix 

oo 

po = | J2 Pfc ' i(T) 1 ^ 1 ® 1 ^ 1 (37) 

where the occupation probabilities of the k and I states are calculated 
assuming, for each atom, a thermal distribution corresponding to temper- 
ature T, as expressed by 



P M (T) = exp 

and 



kn J- 



(38) 



Z=J2 P kAT) (39) 
is the canonical partition function. 



6.1 Gate fidelity 

The most general logical input state has the form 

i 

ix) = c "0\ a >®> ( 4 °) 

a, 0=0 

which is an arbitrary superposition of all two-qubit computational basis 
states. The goal of gate operation is to obtain the ideal output 

l 

\x) = Yl <= a ^\a,0). (41) 

a, (3 = 

This is equivalent to the desired two-qubit transformation Eq. Q: pro- 
vided that i? = 0oo + 0ii — 0oi — 0io, the one can be recovered from the 
other by redefining the logical states via qubit rotations. 

Since in this case the two atoms are described by a density matrix, we 
cannot use the overlap fidelity condition Eq. (| 1 7^ to estimate the perfor- 
mance of the phase gate. We use therefore the minimum fidelity |14| to 
characterize the quality of the phase gate performance 

F = minF(x) 

x 

= min(x | Tr cxt [WS , (p int ® po)S t W t ] | X ) (42) 

x 

Here pint =| x)(x I i s t ne density matrix of the internal state Eq. © and 
S denotes an operator that simmetrizes the atomic state. If we write 

3 

\ x ) = J2cn\n). (43) 

n = 



13 



and assume 

U[\n)®p]tt\n)®U[p], (44) 
the fidelity takes the form 

F = minV \c n \ z \c k \ 2 T nk , (45) 

n,k 

where 

Tnk = e^**-*") P^ n 2iT){n 2 ,n 1 \uiU k \n 1 ,n 2 ). (46) 

The minimum of the fidelity F(x) is evaluated in the Appendix. For the 
ideal case, that is, without anharmonicity but with the exact solutions 
given by (50) , the fidelity at T = is F » 0.99. 

In Tabled we show the values of the fidelity at zero Kelvin with the 
frequencies given in the Table Q which optimizes the fidelity, for different 
values of the anharmonicity. It is important to note that these results are 
obtained considering r = T osc instead of r = 7T oac as in Ref.|5], in order 
to improve the gate time operation. 



A 


F 


1(T 4 


0.99 


lcr 3 


0.96 



Table 2: Fidelity at T = for two values of anharmonicity. 

When the atoms are at finite temperature, the fidelity decreases. If 
we define 7 = exp [— Tiuio/ (fcsT)] and evalute 1 101 neglecting terms of 
the order O (7 s ) the fidelity turns out to be F w 0.97 for A = 10~ 4 
at T ~ 0.5/iK. In Fig. 0we show the behavior of the fidelity with the 
temperature for A = 10 -4 and A = 10 -3 . 

1 



0.8 
0.6 
0.4 
0.2 



"0 °- 1 „ „ <VSl. , 0.3 0.4 

2>c K B T7ih ro ) 

Fi gure 6: Fidelity F against temperature kgT/ (Hujq^X)) for ^"^Rb. Trap parameters have 
the values given in Tabic HI 

It is important to highlight that, since different values of A lead to 
different optimal trap frequencies wo(A), the curves in the Fig.|^]are plot- 
ted as a function of the ratio fcsT/(?iaJo(A)). Thus, e.g., the maximum 
value of the x axis of the A = 10 -4 curve corresponds at a temperature of 
the order of about 1 /iK and for A = 10 -3 is 0.3 (iK, almost an order of 
magnitude smaller. 
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7 Estimating A in a realistic situation 



We show how A is related to the trap's parameters and the properties of 
the atoms used in actual experiments where magnetic traps are used for 
the confinement of neutral atoms. The interaction between the magnetic 
dipole moment of an atom in some hyperfine state \F, tuf) and an external 
magnetic field B is 

flint = -n ■ B. (47) 

In an inhomogeneous magnetic field, if the atomic motion is slow as com- 
pared with the velocity of change of the field vector as seen by the moving 
atom, the interaction only depends on the absolute value of the field: 

Hint = ~^zB = g F m F fJ.BB, (48) 

where \lb is the Bohr magneton and qf is the Lande factor. As in 2 , 
we consider here an atomic mirror like the one realized I10| from a 
solid-state magnetic medium with permanent sinusoidal magnetization 
M = (Mo cos [kMX, 0, 0]) along the x axis. In order to avoid trap losses, 
due to spin flips occuring at magnetic field zeros it is necessary to 
apply a certain external bias field B2 along the y direction. Moreover, 
to obtain a corrugation in the magnetic field modulus, we add a rotating 
external field B\ in the xz plane, at an angle 9 with the surface that can 
be varied at will. In this case the magnetic trapping potential is 

Vm F (x) = g F fJ-Bm F I [B e~ zkM cos 2 (k M x) + f?i cos 0] 2 + 

[B e~ zkM sin 2 (k M x) + B x sin6»] 2 + B^^ . (49) 

where Bo — fXoMo (l — e~ fcMl5 ) /2 and <5 is the tape thickness. The minima 
of Vm F (x) form a periodic pattern above the tape surface, at a height 
zo = fejwr ln(Bo/Si). The spacing between two nearest minima along x 
is just the period of the magnetization 8xm = 2ir/kM- From the series 
expansion of Eq. I|490 along x for 8 = around its minimum, we obtain 
for the anharmonicity parameter 

_ 7T 3 hk M (B x B 2 \ 

A - vm^m \b; + sbJ (50) 

which, for 2nk~^j around a few /tm and for fields B; of a couple hundred 
Gauss, is of the order of 10 -3 . 



8 Conclusions 

In the present paper we have extended the investigations concerning the 
performance of a phase gate, as proposed by Calarco et al. in Ref.|5]. The 
phase gate employs cold trapped neutral atoms and the gate operation 
is obtained with internal state-selective trapping potentials that allow 
collisional interaction between the atoms. 

A correct performance of quantum gates is an essential ingredient of 
a quantum computer. We have therefore relaxed the ideal conditions in 
Ref.|5] in order to check the tolerance of the proposed scheme to exper- 
imental imperfections. We have considered the effects of two possible 
sources of undesired disturbance, i.e., non-perfectly harmonic trapping 
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potentials and temperature. The most crucial parameter is the trap an- 
harmonicity A. By studying the dependence of the gate fidelity on such 
parameter, we have been able to give a prescription to adjust other trap 
parameters as to compensate for this source of infidelity. However, we 
found a critical value for A around 1CP 3 where the fidelity starts to be 
significantly degraded (i.e., well beyond any threshold for fault-tolerant 
quantum computation). This value turns out to be right on the edge 
of what can be presently achieved with magnetic micropotentials (atom 
chips). Thus a further optimization is needed. 

The gate performance could be improved by changing the shape of the 
trapping potential during the gate operation. This requires the simulta- 
neous variation of various physical parameters, A, ui, ujq, and uij_, since the 
fidelity depends on all these quantities. A viable approach to numerically 
search for improved solutions is given by quantum optimal control theory 
|15j . and will be the subject of future investigations. 
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APPENDIX: FIDELITY 



First we see how it is possible to calculate the minimum of the expression 
14511 . Let us define the function 



/3a: 



(51) 



where x~, = |c 7 j 2 with the constraint given by the set of zeros of the 
function 

5(a; 7 ) = ^x„-l. (52) 



Thus we have to solve the linear system of equations 

V£ - AVS = 0. 

that is, 

M a /3Xf3 — A M@a = Tp a + T a f3. 

The minimum of C is then 



(53) 
(54) 



/3a 



a/3 



a/3 



(55) 
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